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ABSTRACT 

Several  numerical  simulations  of  the  transient  flow  of  helium  in  an  expansion  tube  axe 
presented  in  an  effort  to  identify  some  of  the  basic  mechanisms  which  cause  the  noisy  test 
flows  seen  in  experiments.  The  calculations  were  performed  with  an  axisymmetric  Navier- 
Stokes  code  based  on  a  finite-volume  formulation  and  upwinding  techniques.  Although 
laminar  flow  and  ideal  bursting  of  the  diaphragms  was  assumed,  the  simulations  showed 
some  of  the  important  features  seen  in  the  experiments.  In  particular,  the  discontinuity 
in  tube  diameter  at  the  primary  diaphragm  station  introduced  a  transverse  perturbation 
to  the  expanding  driver  gas  and  this  perturbation  was  seen  to  propagate  into  the  test  gas 
under  some  flow  conditions.  The  disturbances  seen  in  the  test  flow  can  be  characterized  as 
either  small-amplitude,  low-frequency  noise  possibly  introduced  during  shock  compression 
or  large-amplitude,  high-frequency  noise  associated  with  the  passage  of  the  reflected  head  of 
the  unsteady  expansion. 


'Research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract 
No.  NASI- 18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science  and 
Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665.  The  author’s  present  address: 
Department  of  Mechanical  Engineering,  The  University  of  Queensland,  St  Lucia,  Qld  4072,  AUSTRALIA. 
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Nomenclature,  Units 

A  :  (x,y)-plane  cell  area 

a  :  local  speed  of  sound,  m/s 

Cp  :  specific  heat  at  constant  pressure,  Jjkg 
Cv  ‘  specific  heat  at  constant  volume,  J/kg 

E  :  total  energy  (internal  +  kinetic),  J/kg 

e  :  specific  internal  energy,  J j  kg 

e  :  unit  vector 

F  :  algebraic  vector  of  a:-component  fluxes 

G  :  algebraic  vector  of  ^/-component  fluxes 

h  :  specific  enthalpy,  J/kg 

k  :  coefficient  of  thermal  conductivity 

n  :  unit  normal  vector 

P  :  pressure,  Pa 

Pr  :  Prandtl  number,  (CppL/k) 

Q  :  algebraic  vector  of  source  terms 

q  :  heat  flux 

R  :  gas  constant,  J/kg/K 

Re  :  Reynolds  number 

r  :  radial  coordinate,  m 

T  :  temperature,  K 

t  :  time,  s 

U  :  algebraic  vector  of  conserved  quantities 
U  :  Riemann  invariant,  m/s 

u  :  x-component  of  velocity,  m/s 

V  :  y-component  of  velocity,  m/s 

X  :  x-coordinate,  m 

y  :  y-coordinate  (radial),  m 

:  parameter  in  Mirels’  boundary  layer  theory 
p  :  density,  kg/rn^ 

7  :  ratio  of  specific  heats 

A  ;  second  coefficient  of  viscosity 

T  :  shear  stress.  Pa 

p  :  coefficient  of  viscosity,  Pa.s 

V  :  frequency, 

n  :  cell  volume,  m^ 

:  volume  per  radian  for  the  axisymmetric  cell 

<  •  >  :  cell-averaged  value 
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Subscripts 

0  :  conditions  immediately  behind  the  shock 

1  ...  20  :  gas  state  as  shown  on  Fig.  1 

e  :  conditions  just  outside  the  boundary  layer 

i  ;  interaction  of  expansion  with  driver-gas/test-gas  interface 

V  :  viscous  contribution 

w,  wall  :  wall  value 

x, y  :  cartesian  components 
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1  Introduction 


The  expansion  tube  is  a  member  of  the  family  of  pulse-type  aerodynamic  facilities  designed 
to  provide  high  enthalpy  test  gas  for  short  durations.  It  has  a  fundamental  advantage 
over  the  reflected-shock  tunnel  because  its  operating  cycle  does  not  involve  the  stagnation 
of  the  heated  test  gas.  Hence,  it  can  provide  a  final  test  flow  with  low  ionization  and 
low  dissociation.  Initial  theoretical  studies  of  the  expansion  tube  [1]  indicated  that  such 
a  machine  could  provide  a  wide  range  of  test  flows  by  simply  altering  the  initial  filling 
pressures.  However,  operational  experience  (see  e.g.  [2,  3])  indicated  that  the  test  flow  was 
often  contaminated  by  large-amplitude  disturbances.  Here,  some  numerical  simulations  of 
an  entire  (but  idealized)  expansion  tube  are  presented  in  order  to  provide  some  insight  into 
the  mechanisms  causing  the  disturbances  in  the  test  flow.  The  specific  facility  configuration 
and  the  flow  conditions  considered  here  are  based  on  the  experimental  perfect-gas  study 
undertaken  by  Shinn  and  Miller  [3]. 

1.1  Facility  Operation 

An  expansion  tube  is  essentially  a  single  tube  divided  into  three  sections  by  “primary”  and 
“secondary”  diaphragms  as  shown  in  lower  part  of  Fig.  1  (which  is  reproduced  from  [3]).  The 
strong  primary  diaphragm  separates  the  driver  tube  and  the  intermediate  (or  shock)  tube 
while  the  light  secondary  diaphragm  separates  the  intermediate  tube  from  the  acceleration 
tube.  The  intermediate  and  acceleration  tubes  have  the  same  diameter.  Initially,  the  driver 
tube  contains  a  “driver  gas”  at  high  pressure  (state  4)  while  the  intermediate  tube  contains 
the  test  gas  (state  1)  and  the  acceleration  tube  contains  a  low  pressure  “acceleration  gas” 
(state  10).  The  aerodynamic  model  to  be  tested  is  located  near  the  downstream  end  of  the 
acceleration  tube. 

The  operation  ot  the  facility,  starting  with  the  rupture  of  the  primary  diaphragm  at  <  =  0, 
is  shown  schematically  in  the  top  part  of  Fig.  1.  The  numbering  of  the  flow  states  is  the  same 
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as  that  introduced  by  Trimpi  [1].  High  pressure  driver  gas  expands  into  the  intermediate 
tube  (state  3)  and  shock-compresses  the  test  gas  (state  2).  Simultaneously  an  expansion  fan 
travels  upstream  into  the  driver  tube.  On  reaching  the  end  of  the  intermediate  tube,  the 
shock  ruptures  the  secondary  diaphragm  (ideally  without  producing  disturbances)  and  the 
leading  portion  of  the  shock- compressed  test  gas  is  allowed  to  expand  into  the  acceleration 
tube  to  reach  the  test  flow  conditions  (state  20).  The  expanded  test  gas  will  have  a  high 
velocity  (typically  6  —  8  kmfs)  and  a  relatively  low  temperature  (several  hundred  degrees 
Kelvin). 

As  shown  in  Fig.  2,  test  time  commences  after  the  passage  of  the  test-gas/acceleration- 
gas  interface  and  ideally  finishes  with  the  arrival  of  the  downstream  end  of  the  unsteady 
expansion  (and  a  corresponding  increase  in  Pitot  pressure).  Not  all  of  the  shock-compressed 
test  gas  is  expanded  to  the  test  flow  state  (20)  and  so  test  times  are  short,  typically  measured 
in  tens  to  hundreds  of  microseconds.  Unfortunately,  ideal  operation  of  the  facility  is  seldom 
realized  as  test  flows  are  typically  very  noisy  and  may  have  peak-to-peak  variations  of  50% 
in  Pitot  pressure.  Note  that,  depending  on  the  length  of  the  tube  sections  and  the  initial  gas 
states,  the  arrival  times  for  the  various  waves  at  the  test-section  may  change.  Of  particular 
importance  is  the  reflected-head  of  the  unsteady  expansion  (of  the  test  gas).  This  wave  is 
generated  when  the  upstream-head  of  the  unsteady  expansion  reaches  the  driver- gas/test- 
gas  interface.  A  partial  reflection  will  occur  and,  for  some  test  conditions,  this  wave  (which 
travels  along  a  u  -f-  a  characteristic)  will  arrive  at  the  end  of  the  acceleration  tube  before  the 
unsteady  expansion  itself. 

1.2  Previous  Studies 

The  cause  of  the  test  flow  disturbances  has  been  the  focus  of  several  studies.  Shinn  and  Miller 
[3]  made  experimental  measurements  in  the  NASA  Langley  facility  using  helium  as  the  test 
gas,  driver  gas  and  acceleration  gas  in  order  to  eliminate  chemical  effects.  They  observed  that 
most  of  the  operating  conditions  resulted  in  test  flows  which  were  very  noisy  but  identified 
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an  optimum  condition  (Pi  =  3.45  kPa,  Pio  =  16  Pa)  with  approximately  300  fis  test  time 
(with  only  moderately  small  disturbances)  before  the  onset  of  large-amplitude  fluctuations 
in  the  Pitot  pressure.  A  large  dip  in  the  Pitot  pressures  was  observed  for  values  of  Pio 
significantly  larger  than  20  Pa  and  this  was  attributed  to  boundary  layer  transition.  They 
also  identified  the  importance  of  the  diaphragm  dynamics  on  the  quality  of  the  test  flow, 
especially  at  low  initial  pressures.  However,  they  did  not  identify  the  mechanism  causing 
the  early  disturbances  in  the  test  gas. 

More  recently,  equivalent  experiments  [4]  using  both  helium  and  argon  as  the  working  gas 
were  performed  at  the  University  of  Queensland  to  check  the  Langley  results  [3]  and  see  if  the 
useful  range  of  test  conditions  could  be  expanded.  In  conjunction  with  these  experiments,  a 
numerical  study  (based  on  quasi-one-dimensional  modelling)  was  undertaken  [5]  to  see  if  the 
test  gas  disturbances  could  be  caused  by  “blobs”  of  driver  gas  being  accelerated  through  the 
expanding  test  gas.  It  was  concluded  that,  even  if  the  blobs  existed,  they  would  not  arrive 
early  enough  to  cause  the  test-flow  disturbances.  However,  the  reflection  of  the  unsteady 
expansion  off  the  driver- gas/test-gas  interface  could  (under  some  conditions)  arrive  early 
enough  to  introduce  the  noise.  Test  time  estimates  based  on  this  mechanism  are  given  in  [6] 
and  are  consistent  with  the  present  numerical  simulations. 

Pauli  [4,  7,  8]  has  pursued  the  possibility  of  acoustic  disturbances  being  generated  near 
the  primary  diaphragm  station  and  then  being  propagated  through  the  unsteady  expansion 
of  the  driver  gas.  A  family  of  transverse  waves  can  exist  and,  when  propagated  through  the 
expansion,  will  undergo  a  frequency  shift  such  that  a  single  dominant  frequency  will  emerge. 
Transmission  of  these  waves  across  the  interface  and  into  the  test  gas  is  assumed  to  be 
determined  by  the  acoustic  impedance  mismatch  at  the  interface.  For  situations  where  the 
sound  speed  in  the  test  gas  is  much  larger  than  that  in  the  expanded  driver  gas,  attenuation 
of  the  waves  in  the  test  gas  should  occur.  Otherwise,  a  significant  fraction  of  the  noise 
will  be  transmitted.  This  acoustic  model  is  seen  to  be  a  reasonable  approximation  to  the 
simulations  presented  here. 
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2  The  Numerical  Method 


The  computations  reported  here  were  performed  with  a  finite-volume  upwind  code  based 
on  the  full  Navier-Stokes  equations.  The  code  is  described  briefly  here  but  further  detail  is 
available  in  [9]. 

2.1  Governing  Equations 


For  an  axisymmetric  flow  (with  y  as  the  radial  coordinate),  the  finite- volume  formulation  of 
the  Navier-Stokes  equations  may  be  expressed  as 


d<U>  1 

dt  q 


+  ^  l(yF  -  yF.)  dy-^,  jj,yG  -  yG.)  dx  =  (?' 


where  <  •  >  indicates  a  volume  average,  Of  is  the  volume  per  radian  and  5  is  a  contour  in 
the  (a:,j/)-plane  ciround  the  volume.  The  £/,  F  and  G  vectors  are 


U  = 


,  yF  =  y 


pu 

pu^  +  P 
pvu 
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puv 

pv^  4-  P 
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Except  for  the  “y”  premultiplying  factor,  these  terms  are  the  same  as  those  in  the  planar 
two-dimensional  situation.  The  viscous  terms  are 
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y^xx 
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Treating  the  viscous  contributions  in  the  form  yr  avoids  any  difficulties  with  the  geometry 
singularity  at  y  =  0.  The  effective  source  term  is 

1  .  V  .( du  dv  v\  ... 

,  where  r„  =  2,r- +  A  +  gj;  +  -  j  .  (8) 

To  augment  these  equations,  the  equation  of  state  for  a  calorically  perfect  gas  is  used 

P  =  p(y-l)e  ,  (9) 

with  7  =  1.667  and  R  =  2077  J/kg/K  for  helium.  The  corresponding  specific  heats  are 
Cv  =  3114  J/kg/K  and  Cp  =  5191  J/kg/K.  The  first  coefficient  of  viscosity  is  evaluated 
as 

H  =  5.023  X  lO-"^  Pa.s  ,  (10) 

and  the  second  coefficient  is  obtained  from  Stokes’  hypothesis  A  =  Also,  a  constant 

Prandtl  number  of  0.67  is  assumed. 

2.2  Numerical  Implementation 

The  flow  domain  in  the  (x,j/)-plane  is  discretized  as  a  structured  mesh  of  quadrilateral  cells 
with  flow  properties  stored  at  the  cell  centres.  At  each  time  step,  the  inviscid-flux  vectors 
(2)  are  evaluated  by  first  applying  a  generalized  MUSCL  interpolation  scheme  [10]  to  obtain 
“left”  and  “right”  states  at  the  midpoints  of  the  cell  interfaces.  A  locally  one-dimensional 
(approximate)  Riemann  solver  is  then  applied  to  obtain  the  interface  flow  properties  during 
the  time  step.  The  spatial  derivatives  used  in  the  viscous  flux  vectors  (3)-(7)  are  obtained 
at  the  cell  vertices  by  applying  the  divergence  theorem.  The  source  term  (8)  is  evaluated  at 
the  cell  centres.  The  line  integrals  in  (1)  are  then  evaluated  using  the  midpoint  rule  and  the 
solution  advanced  in  time. 

Although  a  predictor-corrector  time-stepping  scheme  was  available  in  the  code,  the  solu¬ 
tions  shown  in  the  following  sections  have  used  Euler  time  stepping  in  order  to  reduce  the 


Q'  = 


0 

0 

(P  -  Te0)A/Q' 
0 
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required  cpu  time  for  each  simulation.  Several  simulations  of  the  initial  shock  compression  of 
the  test  gets  were  run  to  investigate  the  effect  of  grid  resolution  and  choice  of  time-stepping 
scheme  [11].  Except  for  the  fine-scale  features,  there  were  few  differences  in  the  solutions  for 
both  Euler  time-stepping  and  predictor- corrector  time-stepping  on  a  (coarse)  1604  x  20  grid. 
The  differences  between  these  coarse-grid  solutions  and  the  solution  on  the  (fine)  2406  x  30 
grid  were  also  small.  The  limit  on  the  number  of  cells  was  set  by  the  memory  available  on 
the  local  workstations  (64  Mbytes)  and  the  desire  to  fit  the  jobs  within  the  8  Mword  queues 
on  a  Cray-YMP  located  at  the  NASA  Langley  Research  Centre,  USA. 
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3  Simulation  Results 


Figure  3  shows  the  flow  domain  considered  in  the  simulations.  The  flow  is  considered  to 
be  axially  symmetric  with  the  flow  domain  extending  radially  from  the  axis  to  the  tube 
wall  and  axially  from  the  closed  end  of  the  driver  tube  to  the  end  of  the  acceleration  tube 
where  it  enters  uhe  test  section  of  the  facility.  Dimensions  are  shown  in  the  figure.  Using 
the  program  in  [11],  the  grid  was  generated  in  three  pieces,  one  for  each  section  of  the 
facility.  The  number  of  cells  and  the  clustering  toward  the  tube  wall  remained  constant 
along  the  length  of  the  domain  and  the  cell  size  remained  constant  in  the  axial  direction. 
The  radial  distance  from  the  outermost  cell-centres  to  the  intermediate-  and  acceleration- 
tube  walls  is  3.87  x  10”"*  m.  Although  the  use  A  an  axisymmetric  geometry  minimizes 
the  computational  resources  required  for  the  simulation,  it  also  precludes  genuinely  three- 
dimensional  disturbances  which  may  constitute  a  significant  portion  of  the  experimentally 
observed  noise. 

While  the  change  in  tube  radius  at  the  primary-diaphragm  station  is  included,  neither 
the  diaphragm  nor  the  square  transition  section  downstream  of  this  point  is  modelled.  In¬ 
stead,  the  simulations  start  with  high-pressure  reservoir  conditions  upstream  of  the  primary 
diaphragm  station  and  intermediate-tube  initial  conditions  downstream  of  the  same  station 
(and  into  the  acceleration  tube).  The  simulations  are  allowed  to  proceed  until  the  primary 
shock  reaches  the  secondary-diaphragm  station  and,  at  that  point  in  time,  the  conditions 
downstream  of  the  secondary  diaphragm  are  reset  to  the  acceleration-tube  initial  conditions. 
Thus,  the  bursting  of  the  diaphragms  is  considered  to  be  ideal  in  the  following  calculations. 
This  is  in  contrast  to  the  experimental  observation  that  the  flow  quality  is  very  sensitive  to 
rupture  behaviour  of  the  diaphragms  [2,  3]  but,  as  will  be  seen  in  the  results,  the  key  features 
governing  the  quality  of  the  test  flow  are  captured. 
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3.1  Initial  State 


All  of  the  simulations  used  the  nominal  driver  tube  conditions  given  in  [3].  These  are 

P4  =  33  MPa,  U4  =  0,  74  =  330  K,  p4  =  48.15  kg/m^,  64  =  1.028  MJ/kg. 

As  discussed  in  [3]  the  nominal  driver-tube  temperature  may  be  lower  than  the  actual  gas 
temperature  at  primary-diaphragm  rupture  where  a  maximum  value  of  T4  =  390  K  was 
estimated.  The  test  gets  and  acceleration  initial  conditions  are 

Tio  =  Ti  =  300  K,  cio  =  Cl  =  0.934  MJ/kg,  Ujo  =  ui  =  0, 

with  pressures  and  densities  as  specified  in  Table  1. 

For  cases  la  and  lb,  the  conditions  were  chosen  to  approximate  the  optimum  conditions 
identified  in  the  experiment  [3]  and  two  different  grids  were  used  to  check  for  sensitivity  to 
grid  resolution.  Cases  either  side  of  this  optimum  condition  were  then  considered.  Case  2  has 
a  higher  intermediate-tube  pressure  where  large-amplitude  noise  was  observed  at  an  earlier 
time  in  the  experiment.  Case  3,  with  a  low  intermediate-tube  pressure,  also  was  observed 
to  have  a  noisy  test  flow  early  in  the  test  time.  Finally,  ceise  4  considers  the  condition  with 
high  acceleration-tube  pressure  where  (for  the  Langley  expansion  tube)  a  large  dip  in  Pitot 
pressure  was  observed  shortly  after  the  arrival  of  test  flow  [3]. 

3.2  Case  1:  Pi  =  3.45  kPa^  Pio  =  16  Pa 

This  condition  was  chosen  to  be  the  starting  point  for  the  study  because  it  was  identified 
as  the  best  operating  condition  in  the  experiment  [3].  The  observed  test  time  w^as  approx¬ 
imately  300  ps  before  large-amplitude  noise  contaminated  the  test  flow.  During  the  test 
time,  moderately  small  disturbances  in  Pitot  pressure  were  observed.  Two  calculations  were 
performed  and  details  of  the  computational  resources  are  shown  in  Table  2.  The  global  time 
step  was  limited  so  that  the  maximum  Courant-Friedrichs-Lewy  (CFL)  number  in  any  cell 
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was  approximately  0.5.  This  maximum  value  occurred  in  the  flow  region  just  behind  the 
shock. 

Ceise  la  used  a  “coarse”  grid  of  1604  cells  axially  by  20  cells  radially  while  case  lb  used 
a  “flne”  grid  of  2406  x  30  cells.  (Refer  to  [11]  for  the  coarse  grid  results.)  The  computed 
results,  including  the  timing  of  events  such  as  the  passage  of  the  shock  and  contact  surfaces, 
were  largely  the  same.  There  were,  however,  some  differences  in  the  details  of  the  contact 
surfaces  and  the  noise  introduced  to  the  test  flow.  In  particular,  the  coarse  grid  attenuated 
some  of  the  noise  and  limited  the  highest  frequency  seen  in  the  simulation. 

Figure  4  summarizes  the  calculation  for  case  lb  as  a  wave  diagram.  The  data  points 
plotted  were  obtained  from  the  individual  plots  of  the  density  field  shown  later.  Note  the 
small  size  of  the  slug  of  compressed  acceleration  gas  (state  20).  The  position  of  the  long 
dashed  line,  denoting  the  downstream-tail  of  the  unsteady  expansion  of  the  test  gas,  was 
estimated  by  its  arrival  at  a;  =  24  m  while  the  short-dashed  line  represents  the  propagation 
of  the  leading  edge  of  the  reflected  expansion  and  the  asterisk  denotes  the  arrival  of  large- 
amplitude  noise  at  the  end  of  the  acceleration  tube  (a;  ~  24  m).  This  reflected  expansion  is 
the  result  of  the  interaction  between  the  upstream-head  of  the  unsteady  expansion  and  the 
test-gas/driver-gas  interface.  The  test  time  (state  5)  for  this  case  starts  after  the  passage  of 
the  test-gas/acceleration-gcis  interface  and  stops  with  the  arrival  of  the  downstream  tail  of 
unsteady  expansion. 

3.2.1  Shock  Compression  of  the  Test  Gas 

Density  contours  for  the  shock  compression  of  the  test  gas  are  shown  in  Fig.  5.  For  each 
frame,  81  contours  have  been  plotted  over  the  range  —6.0  <  logjoP  <  2.0.  Although  the 
shock  heis  remained  essentially  planar,  th’  driver-gas/test-gas  interface  has  become  quite 
distorted  by  <  =  2.0  ms.  Boundary  layer  effects  arc  present  (as  can  be  seen  by  the  contact 
surface  being  left  behind  along  the  tube  wall)  but  they  are  not  well  resolved  in  this  part  of 
the  calculation. 
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Of  special  interest  is  the  train  of  transverse  waves  following  the  interface.  These  transverse 
waves  show  up  clearly  in  the  contours  of  the  radial  velocity  as  shown  in  Fig.  6  where  it 
appears  that  the  disturbances  are  intioduced  at  the  primary  diaphragm  station.  Note  that 
the  axial  scale  is  highly  compressed  and  the  actual  Mach  angles  are  much  shallower  than 
those  that  appear  in  the  figure.  Significant  noise  is  present  in  the  shocked  test  gas. 

Figure  7  shows  the  axial  variation  of  flow  properties  at  f  =  1.0  ms  together  with  the 
estimated  properties  for  a  purely  one-dimensional  situation.  Except  for  some  noise,  the 
levels  are  in  general  agreement.  Although  the  axisymmetric  simulation  exhibits  a  higher 
shock  speed  than  the  one-dimensional  case,  tables  3  and  4  show  that  the  computed  shock 
speed  is  still  lower  than  that  measured  in  the  experiment.  The  large  experimental  value  was 
discussed  in  [3]  and  may  be  partly  due  to  a  low  estimate  for  the  initial  driver-gas  temperature. 
Also,  note  the  variation  in  shock  speed  over  this  phase  of  the  facility  operation  (see  table  4). 

3.2.2  Unsteady  Expansion  of  the  Test  Gas 

The  unsteady  expansion  of  the  test  gas  is  shown  in  Fig.  8.  Viscous  effects  are  more  important 
in  the  acceleration  tube  and  a  very  thick  boundary  layer  is  established  behind  the  shock  and 
test-gas/acceleration-gas  interface.  The  separation  distance  of  the  shock  and  contact  remains 
essentially  constant  from  t  =  2.5  ms  to  t  =  4.0  ms  because  a  balance  has  been  reached 
between  the  mass  flow  of  acceleration  gas  entering  through  the  shock  and  the  mass  flow 
bleeding  into  the  boundary  layer  [12].  The  hot  acceleration  gas  passing  into  the  boundary 
layer  tends  to  accumulate  just  upstream  of  the  test-gas/acceleration-gas  interface  and  form 
a  large  bulge  in  the  boundary  layer  (at  i  ~  23  m,  t  =  4.5  ms).  It  is  not  until  the  last  frame 
(t  =  5.0  ms)  that  any  significant  level  of  noise  is  evident  in  the  expanding  test  gas. 

Figure  9  shows  the  axial  variation  of  flow  properties  at  f  =  4.4  ms.  There  are  very  obvious 
disturbances  in  the  region  of  the  expanded  driver  gas  (state  3)  and,  on  close  inspection  of  the 
pressure  distribution,  some  small-amplitude  disturbances  can  be  seen  near  the  downstream 
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end  of  the  unsteady  expansion  of  the  test  gas  (at  x  22  m).  Pauli’s  acoustic  wave  model  [7] 
suggests  that  the  broadband  disturbances,  introduced  upstream  of  an  unsteady  expansion, 
will  tend  to  focus  to  a  single  frequency  near  the  downstream  tail  of  the  expansion.  For 
the  expanding  driver  gas,  Pauli  estimates  a  “focus”  frequency  vp  =  11.9  kHz  while  the 
perturbations  seen  in  the  present  simulation  have  an  approximate  frequency  v  =  10.2  kHz. 
These  figures  are  in  reasonable  agreement. 

The  very  high  temperature  of  the  shocked  acceleration  gas  enhances  the  diffusion  near 
the  test-gas/acceleration-gas  interface.  Figure  10  shows  details  of  the  shock  and  test- 
gas/acceleration-gas  interface  region.  The  bulge  seen  very  clearly  at  x  ~  23  m,  f  =  4.5  ms 
in  Fig.  8  is  no  longer  as  obvious  in  this  “true-shape”  plot.  Plotting  the  velocity  field  in  a 
frame  of  reference  which  is  stationary  with  respect  to  the  shock  and  contact  surface  clearly 
shows  the  bleeding  of  shocked  acceleration  gas  out  of  region  20  and  into  the  boundary  layer. 
Estimates  of  the  separation  distance  between  the  shock  and  contact  surface  (i.e.  the  length 
of  region  20)  can  be  made  using  the  theory  presented  in  [12].  Using  the  expression 


L20  = 


4^2 


(Pe^Q  \  ^€,0  ^e,0  Pw^O 

PwjO )  ^e,0  Pw,0 


(11) 


together  with  the  parameter  estimates 


/3  =  2.12,  =  Ua,\Q  =  5900  m/s,  Ue.o  =  1200  m/s. 


Pcfi  =  1.0  X  10  kg/wr,  p^jo  =  2.484  x  10  kgjm^,  =  2.012  x  10  ®  Pa.s, 


a  separation  length  of  0.20  m  is  obtained.  This  is  in  reasonable  agreement  with  the  0.22  m 
length  observed  in  Fig.  10  but  it  should  be  noted  that  the  value  of  L20  obtained  from  eqn. 
(11)  is  very  sensitive  to  the  estimated  parameters. 

3.2.3  Histories  of  the  Test  Flow  Properties 

Histories  of  the  centreline  test-flow  properties  near  the  exit  of  the  acceleration  tube  are 
shown  in  Fig.  11.  Features  that  can  be  identified  include  the  passage  of  the  shock  and 
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the  test-gas/acceleration-gcis  interface  at  times  t  =  4.45  ms  and  t  =  4.48  ms  respectively. 
These  events  are  followed  by  a  gradual  variation  in  properties  possibly  due  to  a  constriction 
of  the  effective  tube  diameter  associated  with  the  passage  of  the  bulge  (identified  in  Fig. 
8).  The  downstream  tail  of  the  unsteady  expansion  arrives  at  t  =  4.63  ms  when  average 
pressure  and  density  begin  to  rise  while  Mach  number  and  velocity  decline.  It  is  in  the 
period  4.6  <  t  <  5.0  that  a  “small-amplitude”  fluctuation  is  evident.  The  period  of  this 
perturbation  is  approximately  40/is,  which  corresponds  closely  to  the  period  observed  in  the 
Shinn  and  Miller  [3]  experiment  and  the  more  recent  calibration  measurements  [13]. 

The  fluctuations  increase  in  both  the  amplitude  and  the  frequency  with  the  arrival  of 
the  reflected  head  of  the  unsteady  expansion  at  <  =  5.04  ms.  This  wave  originates  with  the 
rupture  of  the  secondary  diaphragm  and  travels  along  a  U2  ~  ^2  characteristic  until  it  strikes 
the  driver- gas/test-gas  interface.  There  is  both  transmission  into  the  driver  gas  and  partial 
reflection  downstream  along  a  u  -f  a  characteristic  thus  admitting  more  disturbances  into  the 
expanding  test  gas.  The  arrival  time  for  this  wave  (at  the  exit  plane  of  the  acceleration  tube) 
may  be  estimated  as  described  in  [6].  The  process  starts  at  a:'  =  0,  =  0  with  the  rupture 

of  the  secondary  diaphragm.  The  upstream  head  of  the  unsteady  expansion  intersects  the 
driver-gas/test-gcis  interface  at  t'  =  Af„  x’  =  (u2  —  02)  A/,.  Refer  to  tables  5  and  6  for 
particular  values.  The  propagation  time  of  the  wave  through  the  unsteady  expansion  may 
be  obtained  by  integration  along  a  w  -f  a  characteristic  for  which  the  Riemann  invariant  is 
I14| 

2  —  2 

U5  + - -a5  =  U  =  u  + - -a  .  (12) 

7-1  7-1 

Within  the  centred  expansion  we  have 


(13) 


and  this  may  be  substituted  into  equation  (12)  to  obtain 
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u 


7-1 
7  +  1 


7  —  1  <' 


(14) 


The  integration  may  be  terminated  at  the  end  of  the  expansion  (a  =  a^)  or  at  the  end  of  the 
acceleration  tube,  whichever  is  sooner.  If  the  reflected  wave  reaches  the  downstream  end  of 
the  expansion  before  it  reaches  the  end  of  the  tube,  it  continues  to  propagate  into  the  test 
flow  region  (20)  with  speed  us  +  05.  The  time  of  arrival  as  computed  with  this  approach 
closely  matches  the  onset  of  large-amplitude  noise  in  the  exit  flow  histories  and  is  indicated 
by  an  asterisk  in  the  wave  diagram  (Fig.  4). 


A  comparison  of  computed  test  flow  properties  with  the  experimental  measurements  of 
Shinn  and  Miller  [3]  is  shown  in  Fig.  12.  The  values  are  taken  from  the  Pitot  histories  150  fis 
after  the  passage  of  the  shock.  The  overall  level  of  the  computed  wall  pressure  is  a  little  low, 
possibly  because  of  the  low  estimate  for  the  initial  temperature  in  the  driver  tube.  Also, 
the  computed  wall  pressure  is  smoother  than  the  experimental  wall  pressure.  This  may  be 
partly  due  to  numerical  dissipation  and  partly  due  to  noisy  measurement  techniques  in  the 
experiment.  The  arrival  of  the  downstream  tail  of  the  unsteady  expansion  can  be  seen  by 
the  rise  in  pressure  0.2  ms  after  the  passage  of  the  shock  in  both  the  simulation  and  the 
experiment. 


The  other  frames  in  Fig.  12  show  Pitot  pressures.  The  computed  Pitot  pressure  at  the 
centreline  is  again  a  little  lower  than  the  experimental  value  and  is  also  less  noisy  for  the 
first  0.4  ms  after  shock  arrival,  possibly  due  to  excessive  diffusion  in  the  numerical  scheme. 
However,  there  is  good  agreement  in  that  two  types  of  disturbances  are  evident;  a  small- 
amplitude  noise  before  the  arrival  of  the  reflected  head  of  the  unsteady  expansion  and  a 
large-amplitude  noise  after  the  passage  of  the  reflected  expansion.  This  is  one  of  the  main 
results  of  this  study. 

Other  features  evident  in  the  Pitot  histories  include  the  large  dip  at  t  =  5.0  ms,  r  = 
52.5  mm  and  the  consistently  low  values  at  r  =  71  mm.  The  large  dip  is  associated  with 
the  passage  of  the  bulge  in  the  boundary  layer  while  the  r  =  71  mm  position  remains  within 
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the  boundary  layer  after  the  passage  of  the  test-gas/acceleration-gas  interface. 


3.3  Case  2:  Pi  =  20.62  fcPa,  Pio  =  16  Pa 

Having  simulated  the  arrival  of  noise  in  the  test  flow  at  Pi  =  3.45  kPa,  it  was  decided  to 
try  a  higher  pressure  of  Pi  =  20.62  kPa  where  the  experimental  measurements  Indicated 
that  large-amplitude  disturbances  appear  very  early  in  the  test  flow.  Fig.  13  shows  the  wave 
diagram  for  this  case.  Note  that  the  downstream  tail  of  the  unsteady  expansion  arrives  at 
the  tube  exit  shortly  after  the  test  flow  begins  and  that  the  reflected  head  of  the  unsteady 
expansion  is  expected  to  arrive  much  later.  Because  of  the  higher  initial  intermediate-tube 
pressure,  the  shock  has  a  lower  speed  than  in  case  1. 

The  first  frame  of  Fig.  14  shows  the  state  of  the  density  field  just  before  the  rupture  of 
the  secondary  diaphragm.  Transverse  waves  are  still  present  in  the  expanded  driver  gas  but 
they  seem  to  be  weaker.  Also,  there  is  some  noise  evident  in  the  compressed  test  gas. 

The  unsteady  expansion  of  the  test  gas  is  shown  in  sequential  frames  of  Fig.  14.  The 
process  appears  similar  to  that  in  case  1  except  that  there  is  very  little  noise  evident  in  the 
expanded  test  gas  (even  at  t  =  5.0  ms).  This  lack  of  noise  is  confirmed  in  Fig.  15.  Although 
perturbations  have  been  introduced  at  the  primary-diaphragm  station,  the  weaker  expansion 
of  the  driver  gas  has  resulted  in  a  weaker  focusing  of  the  noise  at  the  downstream  end  of  the 
expansion  (when  compared  to  that  in  case  1). 

Fig.  16  compares  the  computed  pressures  at  the  exit  of  the  acceleration  tube  with 
experimental  measurements.  The  arrival  of  the  downstream-tail  of  the  unsteady  expansion 
can  be  clearly  seen  at  <  ~  4.8  ms  while  the  reflected  head  of  the  unsteady  expansion  is  not 
expected  to  arrive  at  the  tube  exit  until  t  ~  5.84  ms  (see  table  6).  The  lack  of  noise  in 
the  laminar  solution  suggests  that  the  disturbances  in  the  experimental  test  flow  are  caused 
by  boundary  layer  transition  or  turbulence  (which  is  not  modelled  in  these  simulations) 
rather  than  propagation  of  noise  through  the  driver-gas/test-gas  interface.  The  omission  of 
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turbulence  is  a  major  limitation  of  this  study  as  the  experimental  flow  is  turbulent  over  most 
of  the  high-pressure  operating  range  of  the  facility. 


The  average  value  of  the  computed  Pitot  profile  (taken  150  ps  after  passage  of  the  shock) 
is  again  slightly  lower  than  the  experimental  value,  possibly  because  of  the  low  initial  value 
for  the  driver  gas  temperature.  The  very  low  Pitot  pressure  values  computed  for  r/r^,  ~  0.7 
are  the  result  of  taking  the  value  of  Pitot  pressure  in  the  hot  acceleration  gas  following 
behind  the  interface. 

3.4  Case  3:  Pi  =  1.74  kPa,  P\q  =  16  Pa 

Given  the  failure  of  the  Pi  =  20.62  kPa  case  to  exhibit  large  disturbances  in  the  simulated 
test  flow,  a  case  with  lower  intermediate-tube  pressure  than  case  1  was  attempted.  Figures 
17  -  19  show  the  results  of  a  simulation  with  Pi  =  1.74  kPa,  Piq  =  16  Pa. 

Table  4  shows  that  the  shock  speed  during  the  initial  compression  of  the  test  gas  is 
slightly  higher  than  for  case  1  but  that  the  final  shock  speed  during  the  unsteady  expansion 
process  is  slightly  lower.  The  shock  compression  of  the  test  gas  appears  to  be  noisier  than  for 
case  1  and  the  centreline  histories  Indicate  that  small-amplitude  disturbances  appear  very 
early  in  the  test  flow.  This  seems  to  be  consistent  with  the  experimental  Pitot  pressure 
results  shown  in  Fig.  19.  The  arrival  of  the  reflected  head  of  the  unsteady  expansion 
{t  ~  5.1  ms)  is  marked  with  a  double  asterisk  in  Fig.  17  and  can  be  seen  as  a  decrease  in 
the  growth  rate  of  Pitot  pressure  in  Fig.  19.  In  the  simulation,  the  arrival  of  this  wave  is 
accompanied  by  larger-amplitude  and  higher-frequency  disturbances  which  are  not  obvious 
in  the  experimental  measurements.  Also,  the  average  value  of  the  computed  Pitot  pressure 
across  the  exit  plane  is  significantly  lower  than  the  experimental  value. 

3.5  Case  4;  P\  =  3.45  kPa,  Piq  =  52.6  Pa 

For  large  values  of  acceleration-tube  filling  pressures,  Shinn  and  Miller  [3]  also  observed  a 
large  dip  in  the  Pitot  pressure  shortly  after  the  start  of  the  test  flow.  The  phenomenon 


17 


was  attributed  to  boundary  layer  transition.  This  particular  case  simulates  such  an  op¬ 
erating  condition  in  order  to  see  if  the  bulge  in  the  boundary  layer  (following  the  test- 
gas/acceleration- gas  contact  surface)  contributed  to  the  large  dip. 


The  results  for  the  laminar-flow  simulation  with  Pi  =  3.45  kPa,  Piq  =  52.6  Pa  are  shown 
in  figures  20  -  22.  The  shock  compression  phase  is  the  same  as  for  case  la  but,  with  the  larger 
mass  of  acceleration  gas,  the  separation  distance  between  the  shock  and  the  contact  surface 
is  increased  and  equilibrium  is  attained  by  ^  =  2.8  ms.  The  histories  of  the  wall  pressure 
shown  in  Fig.  22  indicate  that  small-amplitude  noise  is  present  in  both  the  acceleration  gas 
and  test  gas  up  until  the  arrival  of  the  reflected  head  of  the  unsteady  expansion  at  ^  ~  4.9  ms. 
This  event  is  denoted  on  the  wave  diagram  (Fig.  20)  by  a  double  asterisk  The  arrival  time, 
estimated  with  equations  (12)  -  (14)  is  denoted  by  a  single  asterisk. 

Although  the  large-amplitude  noise  after  the  passage  of  the  reflected  expansion  seems  to 
be  consistent  with  the  experimental  measurements,  there  is  no  evidence  of  the  large  Pitot- 
pressure  dip  in  the  laminar  simulation.  This  is  consistent  with  the  hypothesis  [3]  that  the 
dip  is  caused  by  boundary  layer  transition.  Note  that  this  case  was  computed  on  a  coarse 
mesh  which  may  result  in  the  attenuation  of  the  high-frequency  components  of  the  noise. 
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4  Concluding  Remarks 

From  the  computational  point  of  view,  it  is  possible  to  perform  reliable  simulations  for 
laminar  axisymmetric  flow  of  a  perfect  gas  in  an  expansion  tube  over  a  limited  range  of 
operating  conditions.  Although  viscous  effects  are  adequately  resolved  in  the  acceleration 
tube,  they  were  not  well  resolved  in  either  the  intermediate  tube  or  the  driver  tube.  This 
is  not  considered  to  be  a  major  problem  as  inviscid  multidimensional  flow  features  seem  to 
dominate  the  performance  of  the  driver  and  intermediate  sections  of  the  facility. 

The  numerical  technique  is  capable  of  capturing  strong  shocks  and  contact  surfaces 
with  large  temperature  jumps,  however  the  contact  surfaces  tend  to  be  more  (numerically) 
smeared  than  the  shocks.  The  result  is  that  the  estimates  of  wave  propagation  times  and 
test  flow  properties  appear  to  be  in  reasonable  agreement  with  experimental  measurements. 
However,  the  issue  of  reliably  capturing  the  high-frequency  noise  is  not  fully  settled  and  fur¬ 
ther  computational  experiments  using  alternative  higher-order  schemes  (such  as  ENO  [15]) 
are  required, 

A  laminar  perfect  gas  flow  was  deliberately  selected  to  keep  this  study  manageable,  how¬ 
ever,  there  is  an  immediate  need  to  perform  simulations  with  nonideal  thermal  and  chemical 
effects  with  turbulence  in  three-dimensions.  Diaphragm  dynamics  is  another  important  item 
which  needs  to  be  considered  for  future  study.  In  the  near  future,  it  is  unlikely  thct  all  of 
these  features  could  be  included  but  the  modelling  of  turbulence  and  the  diaphragm  dy¬ 
namics  are  key  items  required  to  simulate  the  expansion  tube  over  a  larger  range  of  useful 
operating  conditions. 

With  respect  to  facility  design  and  operation,  the  simulations  show  some  important 
features.  Firstly,  the  discontinuity  at  the  primary  diaphragm  station  is  sufficient  to  produce 
disturbances  which  may  then  propagate  into  the  test  flow.  Secondly,  the  simulations  clearly 
identify  two  categories  of  noise  in  the  test  flow.  One  is  a  small-amplitude,  low-frequency 
disturbance  which  may  be  seen  before  the  arrival  of  the  reflected  head  of  the  unsteady 
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expansion  while  the  other  is  a  large-amplitude,  high-frequency  noise  associated  with  the 
passage  of  the  reflected  expansion. 

This  distinction  explains  the  experimental  observation  [4,  5]  that  increasing  the  acceleration- 
tube  pressure  Pio  while  keeping  the  intermediate-tube  pressure  P\  constant  results  in  higher 
amplitude  and  higher  frequency  noise  in  the  test  flow.  The  increase  in  Piq  results  in  a  lower 
shock  speed  in  the  acceleration  tube,  while  the  fixed  Pi  means  that  the  arrival  time  of  the 
reflected  expansion  will  be  essentially  constant.  Thus,  at  high  values  of  Pw,  the  reflected 
expansion  is  expected  to  arrive  very  soon  after  the  shock  and  introduce  high-frequency  dis¬ 
turbances  to  the  test  flow. 

To  avoid  the  large-amplitude  noise,  it  should  be  relatively  simple  to  select  operating 
conditions  or  facility  configurations  which  delay  the  arrival  of  the  reflected  expansion.  For 
example,  the  intermediate  tube  could  be  lengthened  while  fixing  the  initial  gas  pressures. 

It  is  not  immediately  obvious  how  to  avoid  the  small-amplitude  disturbances  observed 
before  the  arrival  of  the  unsteady  expansion.  As  can  be  seen  in  Fig.  23  these  disturbances 
do  not  appear  in  the  test  flow  for  all  operating  conditions.  The  acoustic  theory  of  Pauli 
provides  some  guide  to  the  correct  selection  of  good  operating  conditions.  The  basic  results 
of  that  theory  are  that  broad-band  noise  is  focussed  to  a  particular  frequency  across  strong 
expansions  and  that  noise  crossing  the  driver-geis/test-gas  interface  is  attenuated  if  the  speed 
of  sound  in  the  test  gas  is  much  higher  than  the  speed  of  sound  in  the  driver  gas.  Both  of 
these  effects  were  observed  in  the  simulations.  The  stronger  driver-gas  expansions  in  cases 
1,  3  and  4  produced  significant  disturbances  in  the  expanded  driver  gas  while  the  jumps  in 
sound  speed  across  the  driver-gas/ test-gas  interfaces  attenuated  the  noise  admitted  to  the 
expanding  test  gas. 
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Table  1:  Initial  States. 


Case 

CFL 

cpu  —  time 
(hrs) 

memory 
{M  words) 

i  final 

{ms) 

grid 

la 

0.5 

12.5 

3.2 

5.4 

1604  X  20 

lb 

0.5-0.6 

42 

6.8 

5.4 

2406  X  30 

2 

0.5-0.6 

41.6 

6.8 

5.3 

24C6  X  30 

3 

0.5 

4'J 

6.8 

5.4 

2406  X  30 

4 

0.5 

6.5 

3.2 

5.4 

1604  X  20 

Table  2:  Computational  resources. 


Table  3:  Experimentally  measured  shock  speeds  and  pressures  from  Shinn  and  Miller  (1978). 


Case 

{mfs) 

(m/s) 

A 

{kPa) 

Ps.pitot 

{kPa) 

lb 

3745,  3374 

6820,  5968 

0.867 

51.7 

2 

3108 

6160 

1.067 

127.5 

3 

3711 

5645 

0.713 

33.0 

4 

3490 

5570 

0.733 

72.0 

Table  4:  Computed  shock  speeds  and  pressures. 
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Table  5:  Flow  properties  across  the  unsteady  expansion  of  the  test  gas. 


Case 

trxipturc 

(ms) 

Ati 

(ms) 

^tpropagate 

(ms) 

^arrival 

(ms) 

lb 

Dl 

2.96 

5.04 

n 

Ira 

3.44 

5.84 

B 

2.00 

0.95 

2.82 

4.82 

B 

2.08 

0.95 

2.95 

5.03 

Table  6:  Event  times  associated  with  the  reflection  of  the  unsteady  expansion. 
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Figure  1:  Conceptual  wave  diagram  for  the  operation  of  the  NASA  Langley  expansion  tube. 
The  numbers  refer  to  the  states  of  the  gases  as  defined  by  Trimpi  (1962). 
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Figure  2:  Ideal  Pitot  pressure  history  at  the  exit  plane  of  the  acceleration  tube. 
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Figure  3:  Expansion  tube  geometry  as  modelled  in  the  simulations.  Dimensions  are  in 
metres. 
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Figure  4:  Wave  diagram  for  case  1.  Q  =  shock;  A  =  driver-gas/test-gas  interface;  +  = 
test-gzis/acceleration-geis  interface;  long  dash  =  downstream  tail  of  the  unsteady  expansion; 
o  =  upstream  head  of  the  unsteady  expansion;  short  dash  =  reflected  expansion. 


Figure  7:  Flow  properties  along  the  axis  during  shock  compression  (cjise  1,  t  =  1.0  ms).  The 
solid  line  denotes  the  ideal  one-dimensional  solution. 
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Figure  8:  Density  contours  showing  the  unsteady  expansion  of  the  test  gas  for  case  1 
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Figure  10:  Close-up  of  the  flow  near  the  shock  and  contact  surface  at  t  =  4.4  ms,  case  1;  (a) 
density  contours;  (b)  velocity  in  a  frame  of  reference  stationary  with  respect  to  the  shock 
and  contact. 
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Figure  12;  Comparison  of  wall  pressure  and  Pitot  pressure  with  experiment  (case  1,  x 
24  m).  The  experimental  data  is  taken  from  Shinn  and  Miller  (1978). 
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Figure  13:  Wave  diagram  for  case  2.  Q  =  shock;  A  =  driver-gas/test-gas  interface;  +  = 
test-gas/acceleration-gas  interface;  long  dash  =  downstream  tail  of  the  unsteady  expansion; 
o  =  upstream  head  of  the  unsteady  expansion;  short  dash  =  reflected  expansion. 


20.62  kPa,  Pi 


Figure  15:  Properties  along  the  axis  at  <  =  4.6  ms  for  case  2:  (a)  density  contours;  (b)  other 
properties  ais  labelled. 
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Figure  16;  Comparison  of  wall  pressure  and  Pitot  pressure  with  experiment  (case  2,  x 
24  m).  The  experimental  data  is  taken  from  Shinn  and  Miller  (1978). 
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Figure  17:  Wave  diagram  for  case  3.  O  =  shock;  A  =  driver-gas/test-gas  interface;  +  = 
test-gas/ acceleration-gas  interface;  long  dcish  =  downstream  tail  of  the  unsteady  expansion 
o  =  upstream  head  of  the  unsteady  expansion;  short  dash  =  reflected  expansion. 
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Figure  19:  Comparison  of  wall  pressure  and  Pitot  pressure  with  experiment  (case  3, 
21  m).  The  experimental  data  is  taken  from  Shinn  and  Miller  (1978). 


Pi  =  3.45  fcPa,  Pio  =  52.6  Pa 


Figure  20:  Wave  diagram  for  case  4.  O  =  shock;  A  =  driver-gas/test-gas  interface;  +  = 
test-gas/acceleration-gas  interface;  long  dash  =  downstream  tail  of  the  unsteady  expansion; 
o  =  upstream  head  of  the  unsteady  expansion;  short  dash  =  reflected  expansion. 
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